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Abstract-This paper presents the performance of a stand-alone 
self-excited induction generator (SEIG) under balanced/ un- 
balanced excitation with balanced RLC and dynamic load. 
Squirrel cage induction motor has been taken as a dynamic load. 
SEIG is driven with fixed pitch wind energy system. An approach 
based on three-phase induction machine model is employed to 
derive dynamic equations of an isolated SEIG under 
balanced/unbalanced conditions of excitation and balanced static 
and dynamic loads. The SEIG model with balanced/un-balanced 
excitation and balanced load has been simulated using 
MATLAB/SIMULINK. 

Keywords-Dynamic induction motor load; wind turbine; 
balanced loads; balanced/un-balanced excitation capacitor; Self- 
excited induction generator 

Nomenclature 
Yds Vqs Vdr Vqr : Stator and rotor d and q axis voltages 



ids iqs idr iqr : Stator and rotor d and q axis currents 

im : Magnetising current 

Te : Electromagnetic torque 

Lm : Magnetizing inductance 

Vld, Vlq :d and q axis excitation voltage per phase 

ild and ilq :d and q axis load currents per phase 

icd and icq :d and q axis capacitor currents per phase 

ppole : pole pair 

\j/m : Magnetising flux 

Rs and Rr : Stator and rotor resistance 

Xls and Xlr : Stator and rotor leakage reactance 

L : Load imductance 

C : Excitation capacitance 

R : Load resistance per phase 

cor : Rotor angular speed in rad/sec 

p : Differential operator 



I. Introduction 

With increasing concern about the environment and the 
depletion of natural resources such as fossil fuels, much 
research is now focused on obtaining new environmentally 
friendly sources of power. Wind energy is environmentally 
friendly, inexhaustible, safe, and capable of supplying 
substantial amounts of power. Electric power generation using 
non-conventional sources is receiving considerable attention 



throughout the world. The use of induction motor as a generator 
has gained considerable importance to harness energy from 
wind. Self excited induction generator has become popular for 
energy conversion from wind and has advantages of reduced 
cost, size, and ruggedness. Self excited phenomenon in 
induction generator has been well documented in [1-3]. The 
process of voltage buildup in an induction generator is very 
much similar to that of a dc generator and a suitable value of 
residual magnetism should be present in the rotor for the 
process of machine excitation. With wound rotor or squirrel 
cage rotor, the induction generators can be driven by constant 
speed or variable speed drives and can have various 
configurations of different types as: constant- speed constant- 
frequency (CSCF), variable-speed constant-frequency (VSCF), 
and variable-speed variable-frequency (VSVF) [4-5]. 

Tremendous literature is available in the journals and 
conference proceedings on steady state and dynamic operation 
of induction generator under different operating conditions of a 
system. The authors have proposed different models of 
induction generators based on d-q reference theory, impedance 
based approach for model development, admittance based 
models, operational circuit based models, and power equations 
based models [6-11]. Many Authors have reported detailed 
analysis and modeling of induction generator with static as well 
as dynamic model of induction generators. Steady state analysis 
being important both from design and operational point of view 
has been presented by many authors [13-22]. These authors 
proposed steady state analysis technique for analysis of 
induction generator using equivalent circuit of the machine, 
nodal admittance technique, iterative technique etc. to analyze 
the performance of the SEIG. A simple and comprehensive 
approach for steady state analysis was presented [19]. Some 
authors utilized symmetrical components theory as an 
optimization based approached for analysis of SEIG [21-22]. 
Many authors presented transient and dynamic analysis of SEIG 
related to the voltage build up process under self-excitation, 
load perturbation, switching, dynamic load, and unbalanced 
excitation [23-40]. The comprehensive overview of self excited 
induction generator and many related issues has been presented 
nicely [41]. 

Authors have started working in the direction of wind 
turbine driven SEIG in early 90' s in [42-45]. The literature is 
expanded by many researchers in this direction and has been 
reported by many researchers recently on wind driven self 
excited induction generator and dynamic analysis of SEIG 
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[46-49]. State space model of self excited induction generator 
has been presented in [50-53]. 

A. Kishore et al. [51] proposed a generalized state-space 
dynamic modeling of a three phase SEIG has been developed 
using d-q variables in stationary reference frame for transient 
analysis. The state space approach has been found better 
representation of transient operation of SEIG. The d-q axes 
stator-rotor current are the functions of machine parameters. 
The solution of such equation has been obtained assuming all 
the non linear parameters as it is. The modelling of excitation 
system under balanced/un-balanced conditions has been 
developed in terms of d-q. Different constraints such as 
variation of excitation, wind speed and load have been taken 
into account and accordingly the effect on generated voltage 
and current has been analyzed. The effect of excitation 
capacitance on generated voltage has been analyzed. 

In this paper work, an attempt has been made for an analysis 
of fixed pitch based wind driven SEIG and its dynamic behavior 
has been analyzed under no-load, with balanced RLC load as 
well induction motor dynamic load under balanced and 
unbalanced excitation. The simulations have been carried out 
developing model in MATLAB/SIMULINK [54]. The results 
of SEIG under balanced/un-balanced loads, balanced/un- 
balanced excitation are presented in section IV. 

II. Seig Model 

Figure 1 shows the d-q axes equivalent circuit of a (SEIG) 
supplying an inductive load. Equation (18) is a classical matrix 
formulation using d-q axes model is used to represent the 
dynamics of conventional induction machine operating as a 
generator [52]. Using such a matrix representation, one can 
obtain the instantaneous voltages and currents during the self- 
excitation process, as well as during load variations. 




Fig.l d-q axes equivalent circuit of SEIG 



A. Mathematical Model 

The complete dynamic model is represented by the set of 
eight differential equations corresponding to variables ids, iqs , 
idr, iqr, Vld, Vlq as shown in equation (18) is the generalized 
state space representation of a SEIG model. That is in the form 
of classical state-space equation. 

The dynamic model of the three-phase squirrel cage 
induction generator is developed by using stationary d-q axes 



reference frame and the relevant volt-ampere equations are as 
descried as [12]: 

[v]=[Rll]+Lp[l] + co\G\l] (1) 
From which, the current derivative can be expressed as: 

p[ihm[v]-m-o>Mi]) 

p[i]=-[LrMiho>Mi]-[vD (2) 

where [V], [i],[R],[L] and [G] defined below: 

[V]=[V ds V qs V dr V qr ] T (3) 

[i]= [ids iqs idr iqr] (4) 

[R]=diag[R s R s R r R r ] (5) 
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The SEIG operates in the saturation region and its 
magnetizing characteristics are non-linear in nature. 
Magnetizing current should be calculated in every step of 
integration in terms of stator and rotor d-q currents as defined 
by: 

L = ^(i d s+i d rf+(i d s+i d r) 2 (8) 

Magnetizing inductance is calculated from the magnetizing 
characteristics which can be obtained by synchronous speed test 
for the machine under test. This characteristics can be defined 
as: 

L m = 0.1407 +0.0014 i m - 0.0012 i m 2 + 0.00005 ij (9) 
Developed electromagnetic torque of the SEIG is: 

T e =(^polJ^)L m (i qs idr-id S iqr) ( 1 0) 

Torque balance equation is: 

T sh aft=T e + j (2/p) pco r (11) 
The derivative of the rotor speed is: 

/?co r =(p/2)(T shaft -T e )/j (12) 

The residual magnetism in the machine is taken into account 
in simulation process as it is necessarily required for the 
generator to self excite. Initial voltage in the capacitor is 
considered as 2 volts for build-up of voltage for excitation for 
SEIG. 

B. Load and Capac iior Mods!ing+- 1 

Load: Here the modeling of load has been dev eloped in 
tenns of d-q reference frame under balanced conditions. The 
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load currents in terms of their respective voltages have been 
discussed in below: 

RLC-Load Modeling: 

ii d =!((l/L)V ld -(R/L)iid -(1/LC) ji ld ) (13) 

i lq =j((l/L)V lq -(R/L)ii q -(1/LC) ji lq (14) 

The SIMULINK model of RLC load is shown in Fig.2. 



1/LC 



Prod 

a 



1/L 



R/L G> 



►CD 



KD 



Fig.2 SIMULINK model of RLC load 

Capacitor model: 

The capacitor model has been developed based on the 
equations (15), (16), (17) and (18) which represents the self 
excitation Capacitor currents and voltages in d-q axes 
representation as: 
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(18) 
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Where: Kl= Ca + (Cb/2), K2 = Ca + (Cc/2) and K3 = 
+ (K2/Cc) 

This model has been considered for excitation. This model 
is used for both balanced excitation and un-balanced excitation. 
If Ca=Cb=Cc=C then it's called as balanced self-excitation 
otherwise it's called as un-balanced excitation. 

The SIMULINK model of excitation capacitor is shown in 
Fig. 3. This model describes the capacitor model for both 
balanced and unbalanced excitation conditions. 

This is the total block diagram of capacitor model. This has 
been implemented by using equations from (15) to (18). This 
consists two subsystems, these sub systems has been explained 
below from Fig. 4.5(a) to 4.5(b). The developed model in 



SIMULINK for excitation is the proposed new model based on 
the derived equations in section 3. 
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Fig. 3 block diagram of capacitor model 
The following assumptions are made in this analysis: 

• Core and mechanical losses in the machine are 
neglected. 

• All machine parameters, except the magnetizing 
inductance (L m ), are assumed to be constant. 

• The rotor should have sufficient residual magnetism. 

• The three capacitor bank should be of sufficient value 

• Stator voltage in terms of d-q, V ds and V qs should have 
some initial voltage i.e 1 volts. 

The complete model of SEIG with inductive load and 
excitation is given in (19) and with capacitor model is given in 
(20). 



Where: 
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K=l/(L m 2 -L s L r ) 
K 1= C a + (Cb/2), K 2 = C a + (C/2) and K 3 = (K^C) + (K 2 /C c ) 

B. Dynamic Load Model 

The dynamic model of three-phase squirrel cage induction 
motor as dynamic load is similar to that of the induction 
generator but the parameters concerned is related to the motor. 
The volt-current equations of a three-phase induction motor in 
the current derivative can be expressed as: 

pl = [^Wm]ARmVmK[g m Vm]i ( 21 ) 

Where [vM], [iM], [RM], [LM] and [GM] can be defined 
similar manner as defined for SEIG. 

The developed electromagnetic torque of the I. M. is as: 

T = 3Pm L (i i -i i ) (22) 

1 eM ^ ,„ M \ qsM L drM 1 dsM 1 qrM ) v 7 
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(19) 



(20) 



The derivative of electromechanical torque of motor speed can 
be expressed as: 

PO> m =P M ^J M \T L -T eM ) (23) 

Here suffixes d, q refers to D and Q axis (in stationary 
reference frame), s and r refers to stator and rotor, m refers to 
magnetizing component. 

Developed electromagnetic torque of the motor is: 

hs^qr) 

Torque balance equation is: 



T=^L(i i 



qs ar 



(24) 



T shaft —T e — j 



r 2^ 



pco r 



The derivative of the rotor speed from (25) is 
r p\(T shaft -Te) 



pco r = 



(25) 



(26) 



v2; ; 

Here Tshaft is the input torque to SEIG from the turbine, j is 
moment of inertia and p is the number of poles. The SEIG 
operates in the saturation region and its magnetizing 
characteristics are non-liner in nature. So the magnetizing 
current should be calculated in every step of integration in terms 
of stator and rotor D-Q currents as: (reference no): 



(27) 



Magnetizing inductance is calculated from the magnetizing 
characteristics which is obtained by synchronous speed test for 
the machine under test (7.5 kW) and defined as: 

L m = 0.1407 +0.0014 i m - 0.0012 i m 2 + 0.00005 i m 3 (28) 
Three phase generator voltages and currents are obtained 
from D-Q axes components using the relation as: Matrices 
similar relation hold for currents ia, ib and ic. The dynamic load 
model of induction motor has been implemented in SIMULINK 
model using MATLAB. The SIMULINK model is shown in Fig. 
4. 

C. Wind Turbine Model 

Fixed pitch wind turbine model has been utilized to drive 
induction motor. The fixed pitch model of the wind turbine has 
been considered with number of blades as 3, blade radius equal 
to 13m, and the gear ratio is taken as 30 with fixed pitch 

as m Since the power coefficient characteristic of wind 

turbine is a non-linear curve that reflects the aerodynamic 
behavior a wind turbine. The characteristic forms the basis for 

the custom turbine model. The non-linear, dimensionless Cp 
characteristic is represented as [9]. 

— 

C p (A,/3) = C l (-^-C 3 /3-C 4 )e A > + C 6 A ( 29 ) 
A- 
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0.035 



(30) 



x { ;t + o.o8/? /? 3 + i 

where, 

d = 0.5176, C 2 = 116, C 3 = 0.4, C 4 = 5, C 5 = 21, C 6 = 0.0068 
The power coefficient function given by (21), 



-21 



+0.735 



P m =0.5p4(0.5176(— - 9.06> * +0.00681)v w 3 
X 

(31) 

The mechanical torque is given by (22), 



lift — +0.735 

f m =0.5p4(0.5i76( 9.06)e 1 

X 



+ 0.0068A)v w 



<Mv u 



(32) 



The Cp curve is obtained developing the wind turbine model 
in Mat lab SIMULINK. The SIMULINK model of wind turbine 
is shown in Fig. 5. The mechanical power output at air density 
of l.lkpa is obtained as shown in Fig. 6. 
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Fig.4. SIMULINK model of dynamic load 
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Fig. 5 Custom wind turbine model (Air Density: 1.1 kpa) 

Simulink Model of SEIG with Static and Dynamic 
Load 



MATAB SIMULINK is powerful software tool which has 
needed utilized for modeling and simulation of SEIG with RLC 
load and dynamic load under balanced and unbalanced 
excitation. The equations of self-excited induction generator 



(7.5 kW, 4 poles, Rs=Q, Rr=0.77 Q, Xlr=Xls=lrnh.J=0.23 
Kg/m2) described above have been implemented in 
MATLAB/SIMULINK. The equation (19) has been 
implemented in subsystem "Induction Generator" whose 
outputs are currents. The MATLAB SIMULINK model is 
shown in Fig. 7 with subsystems component shown in color. 
Equation (19) shows the eight first order differential equations, 
for which the solutions gives the four currents (stator d-q axis 
currents and rotor d-q axis currents), load currents and capacitor 
voltages. Further these currents are the function of constants viz. 
stator and rotor inductances, resistances, speed, excitation 
capacitance and load impedance. And also variables like 
magnetizing inductance, magnetizing currents and 
electromagnetic torque generated, has been evaluated using 
equations described in the section II. The constraints of non 
linear magnetizing inductance have been taken into account. 
The wind turbine is implemented in subsystem 'wind turbine' as 
shown in SIMULINK model of SEIG in Fig. 7 by using 
equations as described in section II. 

IV. Results and Discussions 

In this work, the results have been determined for SEIG 
(under no-load, with RLC balanced load under balanced and 
constant unbalanced conditions of excitation. In addition to 
RLC load, a dynamic load of induction motor has also been 
considered for study of behavior of SEIG with balanced and 
constant un-balanced excitation. 

The simulations have been obtained using 
MATLAB/SIMULINK model developed for SEIG with RLC 
load as well induction motor load as dynamic load using 
MATLAB 7.0 version. This model has been discussed in the 
previous section. The residual magnetism in the machine is 
taken into account in simulation process as it is necessary 
required for the generator to self excite. Initial voltage in the 
capacitor is considered as 1 volt for build-up of voltage for 
excitation for SEIG. The results of SEIG under balanced and 
dynamic load with balanced/constant un- balanced excitation 
have been determined. Results have been also been obtained 
under no load condition and are discussed in the next section. 
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Fig. 6 Output Mechanical Power of Turbine versus the Turbine Speed 



IJ ES Vol.1 No. 3 2011 PP.151-161www.ijesci.org ©World Academic Publishing 

~155~ 



International lournal of Energy Science 



IIES 



0- 

Vds.Vqs 



wind speed (m/s) 
gen speed 



iqs 




ild 






ids 




L 





capacitor model 



VI d 


Va 




Vb 


VI q 


Vc 


ids 


la 




lb 


iqs 


Ic 


ilq 


ila 




Mb 


ild 


ilc 






R-BALA.NCED 




ilectromannetic torque 



Fig.7. Complete simulink model of wind driven SEIG with load 



A. SEIG Under No-load 

When generator is excited with capacitance value C=l 10 uF 
and rotor speed increased to 310 rad/sec from 50 rad/sec, the 
generated voltage and current attains their steady state value of 
410 volts and 15 amps in 4.6 sec. Cp v/s TSR curve has been 
shown in Fig.8 (graphs 1-4). 

It is observed from the graphs that stator voltage and stator 
currents are balanced and attain their steady state value at 5 sec. 
Rotor speed attains its steady state value in 1 sec. The variation 
of coefficient of power (Cp) vs. tip speed ratio (TSR) for wind 
turbine is shown in Fig. 8 (graph 2) 
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Fig.8 Graphs (1. rotor speed 2. Cp v/s TSR curve 3. stator voltage 4. stator 
currents) 

B. SEIG Operation with RLC Load 

In this case, the behavior of the system connected to RLC 
load is studied. The results have been obtained for two different 
cases. The cases are: 

Case-1: RLC load with balanced excitation 

Case-2: RLC load under constant un-balanced excitation 

Thus, using (3.25), and (3.26), the operation of the SEIG is 
simulated with R=80D, L= 12mH, C= 20 (if under balanced 
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conditions. Generator is excited with C=110 (if and unexcited 
with Ca= 1 10 (if, Cb= 80 of, Cc= 103 uf. 

Case-1: RLC-balanced load with constant balanced excitation 

In Case 1: balanced exaction with balanced RLC-load all 
the currents, i.e. stator currents, capacitor currents and load 
currents are balanced. Voltage is also balanced as shown in 
Fig. 8. It is observed from the Fig. 8 (graphs 1, 2 3, 4) that the 
stator voltage, currents, load currents, and capacitor currents 
attain their steady state value at 3.5 sec. The variation of 
electromagnetic torque with time and magnetizing inductance 
with time of SEIG has also been shown in Fig. 8 (graphs 5 and 
6). The electromagnetic torque is zero initially and then it 
increases exponentially and attains steady values at about 4.2 
sec. It is observed that the magnetizing inductance is constant of 
value 0.14 H up to 3.2 sec. and then reduces to 0.081 H and 
become constant at 4.5 sec. and remains constant. 
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Fig. 8. SEIG at RLC-load with constant excitation (graphs: 1. Stator line 
voltages 2. Stator currents 3. Load currents 4. Capacitor currents. 5. 
Electromagnetic torque 6. Magnetizing inductance) 

Case-2: RLC-balanced load with constant unbalanced excitation 

In Case 2: constant unbalanced excitation and balanced 
RLC load, the load currents are balanced but stator currents and 
capacitor currents are un- balanced as shown in Fig. 9. It is 
observed from the Fig. 9 (graphs 1, 2, 3, 4) that the stator 
voltage, currents, load currents, and capacitor currents attain 
their steady state value at 4.3 sec. The variation of 
electromagnetic torque with time and magnetizing inductance 
with time of SEIG has also been shown in Fig. 9 (graphs 5 and 
6). The electromagnetic torque is zero initially and then it 
increases exponentially and it is observed that Te is having 
oscillations due to unbalanced excitation. It is observed that the 
magnetizing inductance is constant of value 0.14 H up to 3.5 sec. 
and then reduces to 0.087 H and become constant at 5.0 sec. 
and remains constant. 

C. SEIG Operation with Dynamic Load 

In this case, the behavior of the system connected to an 
induction motor load is studied. The results have been obtained 
for two different cases. The cases are: 

Case-1: Dynamic load with balanced excitation 

Case-2: Dynamic load under constant un- balanced 
excitation 

The study of SEIG with induction motor as a dynamic load 
has been studied and implemented using equation as described 
in section II. The operation of the SEIG is simulated with 
induction motor under balanced excitations and constant 
unbalanced excitation. Generator is excited with C=110 (if and 
un-balanced excitation has been taken as with Ca=110 (if, Cb= 
80 (if, Cc= 103 (if. The results obtained with dynamic load are 
shown in Fig. 6. 10 to 6. 1 1 . 

Case-1: Dynamic load with balanced excitation 
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In Case 1 : balanced exaction with balanced dynamic load, 
all the currents, i.e. stator currents, capacitor currents and load 
currents are balanced. Voltage is also balanced as shown in Fig. 
10. It is observed from the Fig. 10 (graphs 1, 2 3, 4) that the 
stator voltage, currents, load currents, and capacitor currents 
attain their steady state value at 5.0 sec. The variation of 
electromagnetic torque of SEIG and motor with time and 
magnetizing inductance with time of SEIG has also been shown 
in Fig. 10 (graphs 5, 6 and 7). The electromagnetic torque of 
SEIG is zero initially and then it increases exponentially and 
attains steady values at about 5.5 sec. The electromagnetic 
torque of induction motor is zero initially and then it decreases 
exponentially and attains steady values at about 5.5 sec. It is 
observed that the magnetizing inductance is constant of value 
0.14 H up to 3.8 sec. and then reduces to 0.095 H and become 
constant at 5.5 sec. and remains constant. 
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Fig.9. SEIG at RLC-load with constant excitation (graphs: 1. Stator line 
voltages 2. Stator currents 3. Load currents 4. Capacitor currents. 5. 
Electromagnetic torque 6. Magnetizing inductance) 
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Fig. 10. SEIG with Dynamic load with constant excitation (graphs: 1. Stator line 
voltages 2. Stator currents 3. Load currents 4. Capacitor currents. 5. 
Electromagnetic torque. 6. Electromagnetic torque of motor 7. Magnetizing 
inductance) 

Case-2: Balanced dynamic load under constant Un-balanced 
excitation 

In Case 2: un-balanced exaction with balanced dynamic load, 
all the currents, i.e. stator currents, capacitor currents are un- 
balanced and load currents are balanced. Voltage is also 
balanced as shown in Fig. 11. It is observed from the Fig. 11 
(graphs 1, 2 3, 4) that the stator voltage, currents, load currents, 
and capacitor currents attain their steady state value at 6.5 sec. 
The variation of electromagnetic torque of SEIG and motor with 
time and magnetizing inductance with time of SEIG has also 
been shown in Fig. 11 (graphs 5, 6 and 7). The electromagnetic 
torque of SEIG is zero initially and then it increases 
exponentially and attains steady values at about 6.5 sec. The 
electromagnetic torque of induction motor is zero initially and 



then it decreases exponentially and attains steady values at 
about 6.5 sec. it is observed that Te is having oscillations due to 
unbalanced excitation. It is observed that the magnetizing 
inductance is constant of value 0.14 H up to 5 sec. and then 
reduces to 0.103 H and become constant at 7 sec. and remains 
constant. 
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Fig.l 1. SEIG at Dynamic load with constant UN excitation (graphs: 1. Stator 
line voltages 2. stator currents 3. Load currents 4. Capacitor currents. 5. 
Electromagnetic torque of SEIG. 6. Electromagnetic torque of motor 7. 
Magnetizing inductance) 

V. Conclusions 

In this paper, the results have been determined for SEIG 
under no-load, balanced RLC load conditions and dynamic load 
using balanced and unbalanced excitation conditions. Based on 
the results for all the cases, the following conclusions can be 
drawn: 

• It is observed that the performance of SEIG under 
balanced RLC load and balanced excitations, the stator voltage, 
stator currents, load currents and capacitor currents are 
balanced. 

• Under constant un-balanced excitation with RLC load, the 
load currents are balanced but stator currents and capacitor 
currents are un- balanced. The electromagnetic torque is 
having oscillations under un-balanced excitation. 

• Under unbalanced load and un-balanced excitation the 
oscillations are more as compared to all other cases. 

• With balanced exaction with balanced dynamic load, all 
the currents, i.e. stator currents, capacitor currents and load 
currents are balanced. They attain their steady values later than 
compared to balanced RLC load. 

• With un-balanced exaction and balanced dynamic load, 
all the currents, i.e. stator currents, capacitor currents are un- 
balanced and load currents are balanced. Voltage is also 
balanced. It is observed that Te is having oscillations due to 
unbalanced excitation. 
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